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Oligomers of heat-shock proteins: Structures that don’t imply function 

William M. JacobsQ Tuomas P.J. Knowles, and Daan Frenkelll 
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Most proteins must remain soluble in the cytosol in order to perform their biological functions. 

To protect against undesired protein aggregation, living cells maintain a population of molecular 
chaperones that ensure the solubility of the proteome. Here we report simulations of a lattice 
model of interacting proteins to understand how low concentrations of passive molecular chaper¬ 
ones, such as small heat-shock proteins, suppress thermodynamic instabilities in protein solutions. 
Given fixed concentrations of chaperones and client proteins, the solubility of the proteome can be 
increased by tuning the chaperone-client binding strength. Surprisingly, we find that the binding 
strength that optimizes solubility while preventing irreversible chaperone binding also promotes the 
formation of weakly bound chaperone oligomers, although the presence of these oligomers does not 
significantly affect the thermodynamic stability of the solution. Such oligomers are commonly ob¬ 
served in experiments on small heat-shock proteins, but their connection to the biological function 
of these chaperones has remained unclear. Our simulations suggest that this clustering may not 
have any essential biological function, but rather emerges as a natural side-effect of optimizing the 
thermodynamic stability of the proteome. 


INTRODUCTION 

Passive molecular chaperones inhibit the aggregation 
of cytosolic proteins and are thus a nearly ubiquitous 
component of living cells (H-Q ■ This class of chaperones 
comprises clusterin, a-crystallins and many other small 
heat-shock proteins (sHSPs), which promote tolerance to 
a wide range of cellular stressors such as elevated tem¬ 
peratures and hazardous nonspecific interactions ii. 
These chaperones cannot by themselves fold or refold 
misassembled proteins and do not require ATP to func¬ 
tion. Instead, passive chaperones associate reversibly 
with aggregation-prone proteins. Even when present 
in sub-stoichiometric ratios with their client proteins, 
sHSPs and similar chaperones are effective at suppress¬ 
ing aggregation and coping with environmental stress [i- 

. Yet the mechanism by which this class of chaperones 
stabilizes the cytosol is not well understood despite sig¬ 
nificant efforts at determining the structural properties 
of these molecules. 

Here we propose that passive chaperones function 
by increasing the overall solubility of the proteome. 
Through this mechanism, passive chaperones reduce the 
fraction of toxic oligomers in solution and suppress the 
nucleation of protein aggregates. It has recently become 
apparent that some sHSPs can also interact with pro¬ 
tein a ggr egates in order to curtail further protein depo¬ 
sition (Md|. These aggregates are often detrimental to 
cellular survival, in part because they can sequester other 
crucial proteins [l2j| . We provide simulation evidence 
that this effect on the proteome solubility is a generic fea¬ 
ture of passive chaperones that associate promiscuously 
and reversibly with their clients. 
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There is substantial experimental evidence that passive 
chaperones interact promiscuously with client proteins in 
chemical equilibrium. Both the rate of client aggregation 
and the fraction of chaperones associated with insolu¬ 
ble proteins are concentration-dependent Si- Further¬ 
more, chaperone binding responds directly to increases 
in the available client binding surfaces, including hy¬ 
drophobic regions of destabilized clients that are only 
transiently exposed [l^. The binding of passive chaper¬ 
ones often modifies the size and structure of amorphous 
aggregates, leading to smaller soluble clusters in which 
the putative chaperone binding sites are associated with 
the hydrophobic interfaces of the client proteins . 

On the basis of these dynamic chaperone-client aggre¬ 
gates, previous studies have suggested that such aggre¬ 
gates might serve as a relatively inert depot of misfolded 
proteins during cellular stress . 

However, client proteins are not the only substrates 
to which passive chaperones bind: these chaperones are 
commonly found in chaperone-only oligomers both in 
vitro and in vivo 0, IT34l^ [2ll - [^ . Recent experi¬ 
ments indicate that these dynamic oligomers are also 
under thermodynamic control [TsI . [l^ . and vary 

with the experimental conditions, such as the tempera¬ 
ture and the ionic strength of the solution [2^, l27l . 
Because this tendency to form oligomers is highly con¬ 
served across the family of sHSPs and similar molec¬ 
ular chaperones, it has long been recognized that dy¬ 
namic fluctuations in the oligomeric state play an im¬ 
portant role in the organization of many passive chaper¬ 
ones 0,111 mill. At present, however, it is unclear 
whether the formation of chaperone oligomers is a key 
functional event. In fact, there is considerable evidence 
to the contrary: experiments have shown that mutations 
and post-translational modifications that alter the ten¬ 
dency of chaperones to form oligomers do not necessar¬ 
ily affect their function [H l3ll - [34l| . These observations 
raise the question of how, if at all, the presence of chap¬ 
erone oligomers contributes to their ability to solubilize 
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aggregation-prone proteins in vivo. Here we show that 
both the function and oligomerization of passive molecu¬ 
lar chaperones can be explained by identifying the opti¬ 
mal conditions for a thermodynamically stable solution of 
chaperones and aggregation-prone proteins. Our results 
suggest that low concentrations of promiscuous chaper¬ 
ones are a generic means of stabilizing a biological mix¬ 
ture with respect to a variety of nonfunctional interac¬ 
tions. 


RESULTS 


Substrate (Client) Chaperone 



Orientationally averaged 
nonspecific interactions: 
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To understand how passive molecular chaperones af¬ 
fect the thermodynamic stability of a protein solution, 
we consider a minimal model of two species in solu¬ 
tion: an aggregation-prone protein and a simple molec¬ 
ular chaperone. Aggregation of the client proteins 
is primarily driven by highly directional interactions. 
These interactions are mediated by ‘patches,’ which 
represent primarily hydrophobic regions that are com¬ 
monly involved in both functional and aberrant protein- 
protein interactions. Chaperone-client recognition is 
also driven by these directional associations between 
chaperone monomers and the exposed patches of client 
monomers. Both the chaperone and client proteins may 
also associate via weak nonspecific interactions, which 
we assume to be averaged over the relative orientations 
of the monomers. These pairwise isotropic interactions 
account for transient associations between proteins in 
a crowded environment [sl, [s^. We do not explicitly 
model the overwhelming majority of proteins that may 
also experience this weak nonspecific interaction but are 
not prone to aggregation via directional interactions, as 
this simplification does not qualitatively affect our anal¬ 
ysis. 


Lattice model of a chaperone—client mixture 

In protein solutions under physiological conditions, the 
interactions between proteins are short-ranged in com¬ 
parison to the size of the monomers, since the high ionic 
strength characteristic of physiological media leads to an 
effective screening of electrostatic interactions [s^]. We 
therefore choose to model protein interactions through 
nearest-neighbor contacts on a three-dimensional lattice, 
where unoccupied lattice sites represent an implicit sol¬ 
vent. Monomers interact if they reside on adjacent lat¬ 
tice sites, and they are free to rotate and to move among 
lattice sites in accordance with the equilibrium Boltz¬ 
mann distribution. We assume that each protein exists 
in a single coarse-grained conformation and that the in¬ 
teractions between proteins are determined by effective 
binding free energies (Figure [T]). This coarse-graining of 
the internal degrees of freedom allows us to capture the 
effects of the intermolecular forces in a reduced set of 
parameters and is particularly suitable for both globular 
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FIG. 1. A minimal model of an associating fluid of passive 
chaperones and aggregation-prone client proteins. Chaper¬ 
one and client monomers interact via nearest-neighbor inter¬ 
actions on a three-dimensional cubic lattice. Orientationally 
averaged nonspecific interactions may be either attractive or 
repulsive. Directional interactions between specific binding 
sites, indicated by blue patches, depend on the relative orien¬ 
tations of the monomers and are always attractive. 


proteins in near-native states and misassembled proteins 
with exposed hydrophobic regions. All monomers on ad¬ 
jacent lattice sites experience an orientationally averaged 
nonspecific interaction, which is assigned a dimensionless 
free energy of —/3e. (Interaction energies are expressed in 
thermal units: I3~^ = /cbT, where fce is the Boltzmann 
constant and T is the absolute temperature.) 

Because aggregation-prone proteins are likely to partic¬ 
ipate in directional protein-protein interactions via mul¬ 
tiple binding sites III IdSl -ldOt. which also promote inter¬ 
actions with sHSPs 1^, we choose a client model 
with three patches that is susceptible to aggregation by 
means of directional interactions alone (Figure [T]). The 
directional interactions between client monomers are as¬ 
signed an attractive free energy of —jdes-s- These inter¬ 
actions are chosen to be strong enough to form insoluble 
client aggregates in the absence of both chaperones and 
additional nonspecific interactions [s^. 

A minimal model of a passive chaperone must be capa¬ 
ble of binding exposed patches on the client monomers. 
Here we assume that the chaperone monomers have a sin¬ 
gle binding site and that the interaction free energy be¬ 
tween chaperone and client patches is —fiecs (Figure [I]). 
While this assumption is clearly a simplification of the 
structure of passive chaperones, which may interact with 
diverse clients via different binding sites, this representa¬ 
tion captures the passivation of interactive client binding 
sites through the burial of hydrophobic surfaces. Most 
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importantly, this representation has the physical features 
that are necessary to capture the qualitative effects of 
passive chaperones on the thermodynamics of a complex 
fluid. Because passive chaperones are known to function 
at low concentrations, we assume that there are always 
fewer chaperones than client binding sites. In what fol¬ 
lows, the relative amounts of the chaperone and client 
monomers in solution are indicated by Xc and Xs , respec¬ 
tively, such that ccc -b Xs = 1. 

We have used Monte Carlo simulations and finite-size- 
scaling techniques to calculate the miscibility limit of 
this model, i.e., the point at which the chaperone-client 
mixture becomes unstable with respect to aggregation 
and/or demixing (see Methods). This miscibility limit 
coincides with a thermodynamic instability, where small, 
spontaneous fluctuations are sufficient to establish long- 
ranged spatial heterogeneity in an initially well-mixed 
solution. In a protein solution, a thermodynamic insta¬ 
bility may have contributions from directional interac¬ 
tions, which cause the polymerization and demixing of 
the strongly interacting species, as well as orientationally 
averaged interactions, which drive the formation of ther¬ 
modynamic phases with differing protein densities [s^. 
Strong directional interactions between the client pro¬ 
teins can thus lead to the formation of disordered aggre¬ 
gates and an accompanying loss of solubility. 


Passive chaperones enhance the thermodynamic 
stability of a protein solution 

As expected, the presence of chaperones inhibits the 
formation of client oligomers by competing for binding 
to patches on the client monomers. However, this passi¬ 
vation of directional interactions is not the only effect of 
chaperone binding: the interactions between chaperones 
and client proteins simultaneously increase the strength 
of the orientationally averaged nonspecific interactions 
that the solution can tolerate while remaining thermo¬ 
dynamically stable. This effect can be seen in Figure [5J 
which shows the miscibility limit, /3e*, at which insolu¬ 
ble aggregates first appear in the solution. When the 
strength of the orientationally averaged nonspecific in¬ 
teractions increases beyond this limit, i.e., /3e > /3e*, the 
solution becomes unstable with respect to small fluc¬ 
tuations in the protein concentrations. Increasing the 
strength of these nonspecific interactions can thus cause 
the solution to become unstable without altering the 
strength of the directional interactions that drive the 
polymerization of the client monomers. Our calculations 
show that passive chaperones dramatically affect the mis¬ 
cibility limit by inhibiting polymerization and solubiliz¬ 
ing transient clusters of client proteins, despite the fact 
that there are far fewer chaperones than there are client 
binding sites. 

In the absence of chaperones, i.e., x^ —t 0, the so¬ 
lution is unstable due to the strong directional interac¬ 
tions between the aggregating client monomers. In this 



FIG. 2. The miscibility limit, /3e*, of a chaperone-client 
mixture depends on the chaperone-client binding strength, 
dcc-s, and the chaperone stoichiometric fraction, Xc- The 
chaperone-chaperone interactions have only a minor effect 
on /3e*: the miscibility limits of solutions with promiscuous 
chaperone interfaces, for which /Icc-c = /3cc-s, are indicated by 
open circles, while the miscibility limits of solutions in which 
chaperone-chaperone interactions are prevented are indicated 
by closed circles. The client-client interaction strength, d^s-s, 
is sufficient to drive the aggregation of clients in the absence 
of chaperones. 

case, I3e* is negative, indicating that a solution of suf¬ 
ficiently concentrated client proteins in a well-screened 
solvent will form insoluble aggregates. It is important to 
note that even when /3ec-s = Oj the chaperones still in¬ 
teract nonspecifically with the client monomers through 
the orientationally averaged interaction f3e. Here we find 
that the addition of such ‘inert’ chaperones has a negligi¬ 
ble effect on the miscibility limit relative to a client-only 
solution. This observation also implies that the major¬ 
ity of cytosolic proteins that are not aggregation-prone 
do not significantly affect the miscibility limit when the 
dominant instability is driven by strong directional inter¬ 
actions. 

Our calculations further indicate that the thermody¬ 
namic forces driving these instabilities are qualitatively 
different in solutions with weakly and strongly binding 
chaperones. In the case of weakly binding chaperones 
(Pcc-s <C /Scs-s), the solution demixes into client-enriched 
and client-depleted phases primarily as a result of di¬ 
rectional interactions. Insoluble client aggregates recruit 
monomers via the formation of directional contacts and 
exchange small oligomers with the coexisting solution. 
With strongly binding chaperones {Pecs ^ P^ss), the so¬ 
lution forms a high-density condensate consisting of both 
chaperones and client proteins bound by nonspecific in¬ 
teractions. Under these conditions, the proteins in both 
the soluble and insoluble phases exist as amorphous clus¬ 
ters that decrease in size as the stoichiometric ratio xdxs 
is increased. The introduction of strongly binding chap- 
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erones, even in low concentrations, significantly increases 
the solution miscibility limit towards the theoretical max¬ 
imum for this model, ~ 0.87. 


Chaperone—chaperone interactions do not 
significantly affect the miscibility limit 

Because a chaperone that interacts with a variety of 
clients is very likely to engage in promiscuous interac¬ 
tions, it is reasonable to assume that chaperones do not 
distinguish among the various hydrophobic surfaces in 
solution. The strength of interactions between chaper¬ 
one binding sites is likely to be similar to the strength of 
interactions between chaperones and clients, and thus it 
is natural to assume that = P^c-s- However, if we 
instead prevent chaperone-chaperone binding by setting 
Pec-c = 0, we find that the effect on the solution misci¬ 
bility limit is negligible (Figure [5]). 

Since the parameter Pec-c directly controls the prob¬ 
ability of chaperone dimerization, our calculations sug¬ 
gest that the formation of chaperone oligomers has a 
very minor effect on chaperone function. Experimentally, 
the relationship between oligomerization and chaperone 
function has been probed by modifying or truncating 
sHSPs [ 23 , . The available experimental evidence 

indicates that alterations to the putative client binding 
sites on sHSPs affect the oligomer equilibria and the func¬ 
tionality of the chaperones independently, in qualitative 
agreement with the predictions of this model. 


A biological fitness function suggests conditions for 
optimized chaperone operation 

Putting these results into context, we now ask, “Is 
there an optimal chaperone-client binding strength for a 
biological mixture?” Figure [5] shows that strongly bind¬ 
ing chaperones are best suited for increasing the mis¬ 
cibility limit. In this case, producing more chaperones 
(or reducing the total concentration of aggregation-prone 
clients) increases /3e* in an approximately linear relation¬ 
ship, allowing an organism to respond effectively to an 
increase in nonspecific interactions. Nevertheless, strong 
promiscuous interactions come at a cost: nearly irre¬ 
versible binding between a chaperone and any available 
association site, including other proteins that are not ex¬ 
plicitly modeled in our simulations, sequesters both in¬ 
terfaces, thereby preventing their participation in further 
functional interactions. The optimal chaperone binding 
strength must balance these competing requirements for 
solution stability and reversible binding. 

Despite the complexity of naturally occurring protein 
solutions, we can predict the optimal chaperone binding 
strength by considering a generic fitness function, which 
quantifies the trade-offs in biological costs and benefits. 
This fitness should be maximized for optimal biologi¬ 
cal function. For the present model, the fitness J-' is 


a function of the miscibility limit, /3e*, as well as two 
biological costs that depend on the variables Pecs and 
Xc- In the absence of any deleterious effects of chap¬ 
erone action, increasing the solubility of the proteome 
must be beneficial, and thus P should be an increasing 
function of /3e*. However, one potential cost of chap¬ 
erone action arises from the sequestration of functional 
proteins (which are not explicitly modeled in our simu¬ 
lations but must be present in a naturally occurring pro¬ 
tein solution) due to the promiscuous binding of chap¬ 
erones. Another potential cost is associated with the 
production of chaperone molecules. These costs imply 
that the fitness function J^[/3e*(/3ec-s,a^c),/3ec-s)a^c] should 
satisfy both dTjdPecs\p^* < 0 and dF/dxc\pc* < 0, re¬ 
spectively, when the miscibility limit Pe* is held constant. 
Taking the total derivative of J- with respect to both Pec-a 
and Xci we find that this fitness function is maximized 
where 


dP£_ 

dPecs 


dT 


dpecs 


dT 

5^ 


and 


dPe* 

dxc 


dT 

dxc 


/3e* 


dT 

'dpi* 


All partial derivatives of T depend on the precise nature 
of the biological system and thus cannot be determined 
precisely. Nevertheless, it is reasonable to assume that 
the cost derivatives of T are approximately constant: at 
low chaperone concentrations, the law of mass action im¬ 
plies that the cost due to reversible, promiscuous bind¬ 
ing is approximately linear in both Xc and Pecs, while 
the total cost associated with the production of chaper¬ 
one molecules is also proportional to their concentration. 
We can therefore interpret the ratios of the derivatives 
in each of the above equations as the importance of each 
cost relative to the benefit of stabilizing the protein so¬ 
lution. Assuming that promiscuous chaperone binding 
and chaperone production are indeed significant biologi¬ 
cal costs, then these equations imply that we should seek 
to optimize the fitness by maximizing the response func¬ 
tions dPe*/dPecs and dpe* jdxc- 

More intuitively, maximizing these response functions 
directs the optimal chaperone design towards the re¬ 
gion of parameter space in which the solution miscibil¬ 
ity limit is most sensitive to small increases in either 
the chaperone-client binding strength or the number of 
chaperone molecules in solution. The first condition, 
dPe*/dpecs, biases the optimal chaperone design away 
from values of Pecs for which the miscibility limit in¬ 
creases asymptotically, thus discriminating against ex¬ 
cessively strong binding between chaperones and clients. 
The second condition, dPe*/dxc, requires that the mis¬ 
cibility limit be sensitive to changes in the chaperone 
stoichiometric fraction. 

Our calculations show that it is indeed possible to sat¬ 
isfy both conditions simultaneously. In Figures [3^ and 
[3b, we plot the calculated response functions dPe*/dPecs 
and dPe*/dXc, in dimensionless units, as functions of 
the chaperone stoichiometric fraction and the chaperone 












FIG. 3. The optimized design window for a passive molecular chaperone coincides with the conditions under which chaperone 
oligomerization is most probable, (a) The response in the solution miscibility limit to an increase in the chaperone stoichiometric 
fraction. The approximately linear response regime for strong binding chaperones is indicated in red. (b) The response in the 
solution miscibility limit to an increase in the chaperone-client binding strength, (c) The probability that a chaperone binding 
interface is bound to another chaperone monomer. For a given stoichiometric ratio of chaperones and clients, this probability 
is greatest in approximately the same design window in which both dfit* /dxc and dPe* /dfitc-s are simultaneously maximized. 


binding strength. We identify a ‘design window’ for op¬ 
timal chaperone operation by finding the approximate 
range of chaperone binding strengths over which both re¬ 
sponse functions are maximized given a fixed chaperone 
stoichiometric fraction. The region of parameter space in 
which both response functions can be maximized is rela¬ 
tively narrow, suggesting that optimized passive chaper¬ 
ones should have tightly constrained binding strengths. 
We further find that the optimal range of chaperone bind¬ 
ing strengths is only weakly dependent on the chaperone 
stoichiometric fraction. This observation implies that 
chaperones with fixed binding strengths can operate close 
to optimality over a wide range of sub-stoichiometric con¬ 
centrations. We also note that these protein-protein 
interaction free energies are in the physical range of a 
few /cbT. The optimal chaperone binding strength is 
generally weaker than the client-client interactions, in¬ 
dicating that the chaperones need not out-compete the 
aggregation-prone clients for association with exposed 
binding interfaces. 


The optimal chaperone—client binding strength 
promotes chaperone oligomerization 

Remarkably, our simulations reveal that the probabil¬ 
ity of finding chaperone oligomers is also highest in the re¬ 
gion of parameter space where the optimal design condi¬ 
tions for chaperone activity are satisfied. In Figure[3]:, we 
plot the probability of chaperone-chaperone binding at 
the miscibility limit, assuming that Pec-c = P^c-s- We find 
that this probability is maximal in the window of opti¬ 


mal chaperone binding strength over the complete range 
of simulated chaperone stoichiometric fractions. Under 
these conditions, a significant fraction of the chaperone 
binding sites are not associated with the aggregation- 
prone interfaces on the client proteins, but are rather 
buried in chaperone-only oligomers. This fraction may 
be even higher in the miscible fluid or in the presence of 
client proteins that are less prone to aggregation due to 
weaker directional interactions. 

These calculations provide further evidence that the 
assembly of chaperone oligomers does not play a func¬ 
tional role. Although the choice of Pcc-c affects the mag¬ 
nitude of the effect shown in Figure [3J:, we emphasize 
that simply allowing chaperone-chaperone binding does 
not imply that chaperone-only oligomers will be observed 
at the miscibility limit: there is a large region of pa¬ 
rameter space over which this probability is very small. 
Furthermore, the results presented in Figure [3] are quali¬ 
tatively unchanged for all reasonable choices of Pec-c, i.e., 
0 < Pcc-c ^ Psc-s- Our simulations thus indicate that the 
ability to assemble chaperone oligomers affects neither 
the anti-aggregation function of the chaperones nor their 
adherence to the proposed design constraints. 


DISCUSSION 

We have shown that a simple model of chaperone- 
client mixtures reveals two generic and unexpected fea¬ 
tures of passive molecular chaperones. First, chaperone- 
chaperone interactions only marginally affect the stabil¬ 
ity of a protein solution in which strong directional inter- 
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actions drive the aggregation of client proteins. Second, 
promiscuous passive chaperones tend to assemble chap¬ 
erone oligomers under conditions where the chaperone- 
client binding strength balances the requirement for pro- 
teome stability with the need to avoid irreversible bind¬ 
ing. Taken together, these results suggest that the assem¬ 
bly of oligomers of passive molecular chaperones is not 
an essential functional event for stabilizing a protein so¬ 
lution. Instead, this behavior emerges as a side-effect of 
operating under thermodynamically optimal conditions. 

To arrive at this conclusion, we have proposed that 
passive chaperones perform their anti-aggregation func¬ 
tion by increasing the miscibility limit of a protein so¬ 
lution. Through this mechanism, passive chaperones in¬ 
hibit the sequestration of functional proteins and increase 
the thermodynamic stability of a biological mixture with 
respect to random nonspecific interactions. Our simu¬ 
lations demonstrate that this mechanism is physically 
plausible even when the aggregation-prone client proteins 
greatly outnumber the chaperones. We emphasize that 
only the ratio of chaperone molecules to client binding 
interfaces, not the total concentration of chaperones in 
solution, is relevant for chaperone function. In all cases 
considered here, the stoichiometric fraction of chaperones 
is much lower than Xc = 0.75, the fraction that would 
be required to passivate all binding sites on the three- 
patch client monomers in the solution. The fact that 
the chaperones that fulfill this anti-aggregation function 
are highly conserved in both lower and higher organisms 
suggests that there is a strong evolutionary pressure to 
perform this role in an optimized fashion. Our calcula¬ 
tions indicate that the range of suitable chaperone bind¬ 
ing strengths is indeed narrow and that the principles 
for an optimal design emerge from thermodynamic argu¬ 
ments. 

The generality of the present model suggests that the 
assembly of chaperone-only oligomers would not be af¬ 
fected by introducing additional detail in an off-lattice 
model. Such an extension, however, would allow a much 
wider variety of chaperone oligomers to be observed. The 
significant coarse-graining involved in the development of 
the present model and the high symmetry imposed by the 
lattice do not permit the reproduction of many structural 
features or the precise oligomeric distributions of specific 
passive chaperones. For instance, all three domains of 
aB-crystallin are believed to be involved in the assembly 
of higher-order oligomers [^ . while the chaperones in 
the present model may only form dimers through client¬ 
binding interfaces. Nevertheless, such detailed molecular 
interactions are unlikely to affect the physical mechanism 
by which passive chaperones suppress aggregation. 

Most importantly, the simplicity of this model allows 
us to make generic predictions about the thermodynam¬ 
ics of passive molecular chaperones. Regardless of the 
molecular-level details, the critical behavior of a fluid 
with short-ranged interactions falls within the Ising uni¬ 
versality class, which is also known to describe phase 
separation in globular protein solutions In the 


vicinity of the miscibility limit, fluctuations in both the 
protein density and the intermolecular contacts within 
aggregates are significant, and a broad distribution of 
cluster sizes is observed at equilibrium. Our proposed 
mechanism therefore supports the assertion that sub¬ 
unit exchange is essential for the function of sHSPs and 
related chaperones i a M m. m. Even if the ag¬ 
gregates are not fully equilibrated due to slow kinetics, 
the large concentration fluctuations in the vicinity of a 
metastable critical point are likely to enhance the forma¬ 
tion of gel-like a ggreg ates [50j | or the nucleation of aggre¬ 
gated phases |5lL l52j ]. For example, recent simulations 
have shown that clustering through nonspecific interac¬ 
tions plays an important role in the kinetics of amyloid 
fibril nucleation 15311. 


CONCLUSIONS 

We have presented a minimal model of a mixture of 
passive molecular chaperones and aggregation-prone pro¬ 
teins. By calculating the limit of thermodynamic sta¬ 
bility in this model protein solution, we have shown 
how passive chaperones that are expressed in sub- 
stoichiometric ratios with their clients can substantially 
suppress aggregation. We have further argued that the 
biological costs associated with chaperone production 
and promiscuous, irreversible binding significantly con¬ 
strain the optimal design of an effective passive chaper¬ 
one. We find that if passive chaperones interact promis¬ 
cuously with exposed hydrophobic surfaces, then the as¬ 
sembly of chaperone oligomers emerges as a nonfunc¬ 
tional side-effect of this thermodynamically optimal de¬ 
sign. Because of the generality of the model, these con¬ 
clusions are relevant to a broad class of molecular chap¬ 
erones. Fully atomistic simulations could provide further 
information on the parameters governing the interac¬ 
tion strengths between chaperones and their aggregation- 
prone targets as well as between the passive chaperones 
themselves. Such simulations could therefore provide a 
means of transferring the general thermodynamic prin¬ 
ciples uncovered by the coarse-grained simulations pre¬ 
sented here to detailed models of specific chaperone- 
client mixtures. 


METHODS 

In the lattice model considered here, the limit of ther¬ 
modynamic stability of a well-mixed solution is encoun¬ 
tered at the critical surface for phase separation. In what 
follows, we describe the Monte Carlo simulations and 
finite-size-scaling theory used to calculate points on this 
critical surface. Our approach is a generalization of the 
computational strategy described in detail in Ref. HI. 

In general, the critical surface of a multicomponent 
mixture has dimension d — 2, where d is the total num¬ 
ber of independent thermodynamic fields . The inde- 
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pendent thermodynamic fields in the present model are 
the dimensionless chemical potentials of both the chap¬ 
erones and the clients, /3^c and /3/is, respectively, as well 
as the dimensionless interaction energies: /3e, Pcg-s, p£c-s 
and /3ec-c- The relevant critical surface in this model is 
thus a 4-dimensional surface. 

We perform biased grand-canonical Monte Carlo simu¬ 
lations, as described in Ref. IH, to collect statistically in¬ 
dependent lattice configurations near the critical surface. 
We use a L X L X L cubic lattice with periodic boundary 
conditions and set L = 12 so that all simulations are car¬ 
ried out in the scaling regime. We then apply th e finite- 
size-scaling theory of Wilding and Bruce [5^ to solve 
self-consistently for the critical order parameter, A4, and 
the critical orientationally averaged nonspecific energy, 
/3e*, at fixed values of /3es-s, /3ec-s, /3ec-c and Xc- In order 
to determine each critical point plotted in Figure [2l we 
approximate the marginal probability distribution p{A4) 
from the grand-canonical samples and then tune this dis¬ 
tribution in order to match the known distribution of the 
critical ordering operator in the three-dimensional Ising 
universality class, pm- This computational procedure is 
described below. 

In a two-solute solution, with two independent dimen¬ 
sionless chemical potentials Ppc and /3/Xs, the critical or¬ 
der parameter must account for fluctuations in the num¬ 
ber densities of both the client and chaperone monomers, 
Ps and Pc, respectively, as well as fluctuations in the in¬ 
ternal energy density, u. The critical fluctuations in the 
number densities can be described by the vector u, which 
indicates the difference in compositions of the two incip¬ 
ient phases [s^. We therefore define Ai to be the linear 
combination 


M = VsPs A VcPc - su, 


( 1 ) 


where both i> and the field-mixing parameter s must be 
determined self-consistently. The grand-canonical distri¬ 
bution of A4 is constructed from the simulation data ac¬ 
cording to 


pj^k = Ay^ < {ps,Pc,u)y-M < 6Mk+i'^ 


( 2 ) 

where the index v runs over all independent samples and 
!{•} is the indicator function. Each sample has a sta¬ 
tistical weight Wy in the grand-canonical ensemble that 
depends on the values of the thermodynamic fields [s^ . 
The system-dependent scaling constant A must be deter¬ 
mined self-consistently. The bin size is chosen such that 
{SMk+i - SMk) = where SM = A{M - M*) and 
A4* is the ensemble-averaged mean value of A4. 


We then construct a x^-function that seeks to minimize 
the difference between the observed distribution of A4 
and the universal distribution, pm, while obeying the 


imposed composition constraint: 


k 


[pi^k iPf) - Pm {SMk/A) 


+ 


( 3 ) 


{(p^{|3f))/Y.Je{s,c}{pJ{|Sf))-Xi) 

iG{s,c} 


where fif = {fie, Pe^-s, Pecs, Pec-c, fips, Ppc) and the in¬ 
dex k runs over all bins. In the second term, {pp in¬ 
dicates the ensemble-averaged number density of compo¬ 
nent i. We estimate the error in the sampled distribution 
of M to be 

2 {^2v^v^k,v) ~ (X/ii /^samples , , 

^k ; v^l 

Hv Wv 


where \k,v is the indicator function written out explicitly 
in Eq. m, and we estimate the error in the observed 
composition at the critical point to be 



E (^Pj^Pk) 



where {5pj5pk) = (pjPk) - {Pj){Pk), P =J2je{s,c}Pj and 
Sij is the Kronecker delta. 

Finally, we calculate the probability of chaperone 
dimerization, (pc-c), directly from the simulation data 
according to the definition 


(6) 

where rice is the number of chaperone-chaperone patch 
contacts and Nc is the total number of chaperone 
monomers on the lattice. In this definition, (•)* indicates 
a grand-canonical average obtained at the critical point 
with the specified chemical potentials and directional in¬ 
teraction energies. 
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